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The recent WMAP data have confirmed that exotic darlc matter together with the vacuum energy 
(cosmological constant) dominate in the fiat Universe. Thus the direct dark matter search, consisting 
of detecting the recoihng rmcleus, is central to particle physics and cosmology. Modern particle 
theories naturally provide viable cold dark matter candidates with masses in the GeV-TeV region. 
Supersymmetry provides the lightest supersymmetric particle (LSP), theories in extra dimensions 
the lightest Kaluza-Klein particle (LKP) etc. In such theories the expected rates are much lower 
than the present experimental goals. So one should exploit characteristic signatures of the reaction, 
such as the modulation effect and, in directional experiments, the correlation of the event rates with 
the sun's motion. In standard non directional experiments the modulation is small, less than two 
per cent and the location of the maximum depends on the unknown particle's mass. In directional 
experiments, in addition to the forward-backward asymmetry due to the sun's motion, one expects a 
larger modulation, which depends on the direction of observation. We study such effects both in the 
case of a light and a heavy target. Furthermore, since it now appears that the planned experiments 
will be partly directional, in the sense that they can only detect the line of the recoiling nucleus, 
but not the sense of direction on it, we study which of the above mentioned interesting features, if 
any, will persist in these less ambitious experiments. 

PACS numbers: 95.35.+d, 12.60.Jv 



INTRODUCTION 

The combined MAXIMA-1 BOOMERANG 0, DASI Q and COBE/DMR Cosmic Microwave 
Background (CMB) observations 0] imply that the Universe is flat [H] and that most of the matter in 
the Universe is Dark i.e. exotic. These results have been confirmed and improved by the recent 
WMAP data [7]. Combining the the data of these quite precise experiments, crudely speaking, one 
finds: 

= 0.05, flcDM = 0.25, flA = 0.70 

Since the non exotic component cannot exceed 40% of the CDM [8], there is room for the exotic 
WIMP's (Weakly Interacting Massive Particles). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the 
observed rotational curves, it is essential to directly detect '^-1^ such matter. Until dark matter is 
actually detected, we shall not be able to exclude the possibility that the rotation curves result from 
a modification of the laws of nature as we currently view them. This makes it imperative that we 
invest a maximum effort in attempting to detect dark matter whenever it is possible. Furthermore 
such a direct detection will also unravel the nature of the constituents of dark matter. The possibility 
of such detection, however, depends on the nature of the dark matter constituents. 

Supersymmetry naturally provides candidates for the dark matter constituents [^-[lH. In the 
most favored scenario of supersymmetry the LSP (Lightest Supersymmetric Particle) can be simply 
described as a Majorana fermion, a linear combination of the neutral components of the gauginos 
and higgsinos In most calculations the neutralino is assumed to be primarily a gaugino, 

usually a bino. 

Since the WIMP's are expected to be very massive, m^^jMP > 30GeV, and extremely non 
relativistic, with average kinetic energy -< T 50KeV{m]vjMP /WOGeV), they are not likely 
to excite the nucleus. So they can be directly detected mainly via the recoiling of a nucleus (A,Z) in 
elastic scattering. The event rate for such a process can be computed from the following ingredients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained in the framework 
of the prevailing particle theory. For supersymmetry this is achieved as described , e.g., in 
Refs 



2. A well defined procedure for transforming the amplitude obtained using the previous effective 
Lagrangian from the quark to the nucleon level, i.e. a quark model for the nucleon. This step 
in SUSY models is not trivial, since the obtained results depend crucially on the content of 
the nucleon in quarks other than u and d. 

3. Knowledge of the relevant nuclear matrix elements — obtained with as reliable as 
possible many body nuclear wave functions. Fortunately in the case of the scalar coupling, 
which is viewed as the most important, the situation is a bit simpler, since then one needs 
only the nuclear form factor. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. Since the essential 
input here comes from the rotational curves, dark matter candidates other than the LSP 
(ncutralino) are also characterized by similar parameters. 

In the standard nuclear recoil experiments one has the problem that the reaction of interest does 
not have a characteristic feature to distinguish it from the background. So for the expected low 
counting rates the background is a formidable problem. Some special features of the LSP-nuclear 
interaction can be exploited to reduce the background problems. Such are: 

• The modulation effect. 

This yields e periodic signal due to the motion of the earth around the sun. Unfortunately 
this effect is small, < 2% for most targets. Furthermore it is inevitable to have backgrounds 
with a seasonal variation. 

• Transitions to excited states. 

In this case one need not measure nuclear recoils, but the de-excitation 7 rays. This can 
happen only in vary special cases since the average WIMP energy is too low to excite the 
nucleus. It has, however, been found that in the special case of the target ^^^I such a process 
is feasible [l^ with branching ratios around 5%. 



Detection of electrons produced during the WIMP-nucleus collision. 

It turns out, however, that this production peaks at very low energies. So only gaseous TPC 
detectors can reach the desired level of lOOeV. In such a case the number of electrons detected 



may exceed the number of recoils for a target with high Z [17| , [1 



• Detection of hard X-rays produced when the inner shell holes are filled. 
It has been found [3] that in the previous mechanism inner shell electrons can be ejected. 
These holes can be filled by the Auger process or X-ray emission. For a target like Xe these 
X-rays are in the 30keV region with the rate of about 0.1 per recoil for a WIMP mass of 100 
GeV. 

In the present paper we will focus on the characteristic signatures of the WIMP nucleus interaction, 
which will manifest themselves in directional recoil exp eriments, i.e. experiments in which the 



direction of the recoiling nucleus is observed [20| , [2]| , [22] , [23| , [2J| , [25[ , [26^ We will concentrate 



on the standard Maxwell-Boltzmann (M-B) distribution for the WIMPs of our galaxy and we will not 
be concerned with other non thermal distributions, even though they may yield stronger directional 
signals. Among those one should mention the late infall of dark matter into the galaxy, i.e caustic 
rings [^, 29, 30, 31, 32], dark matter orbiting the Sun [i^] and Sagittarius dark matter [s^. 

We will present our results in such a fashion that they do not depend on the specific properties 
of the dark matter candidate, except that the candidate is cold and massive, m^/i/iMP ^ 10 GeV. 
So the only parameters which count is the reduced mass, the nuclear form factor and the velocity 
distribution. So our results apply to all WIMPs. In a previous paper we have found that the 
observed rate is correlated with the direction of the sun's motion 2^, 13] ■ On top of this one will 



observe a time dependent variation of the rate due to the motion of the earth. Those features 
cannot be masked by any known background. Unfortunately, however, even the most ambitious of 
the planned experiments are not expected soon to distinguish the two possible senses along the line 
of nuclear recoil [s^]- Such experiments cannot, e.g., measure the backward-forward asymmetry. 



on this occasion we will extend our previous directional calculations [2g|- 27| and explore, which 
characteristics, if any, of the previous calculation persist, if the rates for both senses of the nuclear 
recoil are summed up. 



RATES 



Before computing the event rates for WIMP nucleus scattering we should discuss the kinematics. 
In the case of the WIMP-nucleus collision we find that the momentum transfer to the nucleus is 
given by 

q — 2/irU cos 6*, (1) 

where 9 is the angle between the WIMP velocity and the momentum of the outgoing nucleus and 
is the reduced mass of the system. Instead of the angle 9 one introduces the energy Q transferred 

2 

to the nucleus, Q = 2 Am (^™p is the nuclear mass). Thus 

2 sin 9 cos 9d9 = --r^^^dQ 
Furthermore for a given energy transfer the velocity v is constrained to be 



^ L- ^min 1 ^mm — A/ o ^ ' 

V ^ fir 

We will find it it convenient to introduce instead of the energy transfer the dimensionless quantity u 

u = liqb)^ ^ AmpQb'' ^ u ^ , Qo = X 10^ A-^/^ KeV (3) 

where h is the nuclear (harmonic oscillator) size parameter. 

It is clear that for a given energy transfer the velocity is restricted from below. We have already 
mentioned that the velocity is bounded from above by the escape velocity. We thus get 



aVu < y < nyesc , a 

with n > 1 (see below). 



V2^irbvQ\ (4) 



d u 

2 sin9 cos 9d9 — ^-dj/, y = — (5) 

y vo 

The differential (non directional) rate with respect to the energy transfer u can be written as: 

dUundir = -. da[u,v)\v\ (6) 

AniN 

Where p{Q) — 0.3GeV/cm^ is the LSP density in our vicinity, m is the detector mass, is the 
WIMP mass and da{u, v) the nucleus WIMP cross section. 

The corresponding directional differential rate, i.e. when only recoiling nuclei with non zero velocity 
in the direction e are observed, is given by : 

p(0) Tn 1 

dRdi,- = \v\v.e Q{v.e) — da{u,v) (7) 

ArriN ^vr 

^/u ... . / 1 ,a; > 



The LSP is characterized by a velocity distribution. For a given velocity distribution with 
respect to the center of the galaxy, One can find the velocity distribution in the lab frame f{v,VE) 
by writing 

Vq is the sun's velocity (around the center of the galaxy), which coincides with the parameter 
of the Maxwellian distribution, and Vi the Earth's velocity (around the sun). The velocity of the 
earth is given by 

Ve — vqz + Wi( sina x — cosa cosjy + cosa sinj z ) (8) 

In the above formula z is in the direction of the sun's motion, x is in the radial direction out of the 
galaxy, y is perpendicular in the plane of the galaxy {y — z x x) and 7 « 7r/6 is the inclination of 
the axis of the ecliptic with respect to the plane of the galaxy, a is the phase of the Earth in its 
motion around the sun (a = around June 2nd). 

The above expressions for the rates must be folded with the LSP velocity distribution. In the 
present work we will assume that the velocity distribution is Maxwell-Boltzmann (M-B) with a 
characteristic velocity vq and an upper cut off equal to the escape velocity, Vesc = 2.84uo. For 
comparison we will also consider M-B distributions with larger characteristic velocity nvo and cut 
off velocity 2.84nt;o, with n > 1. This situation arises, if one considers the coupling of dark matter 
to dark energy via a scalar field. Then the gravitational interaction of matter is not affected. The 
interaction involving dark matter is increased. Via the virial theorem this results to an isothermal 
M-B distribution with higher temperature. We will not elaborate further on this point, but we refer 



the interested reader to the literature. [35|, [3 
We will distinguish two possibilities: 



The direction of the recoiUng nucleus is not observed. 

Even though our main interest is in the directional rate for orientation purposes we will summarize 
the main points entering the standard (non directional) searches. The non-directional differential 
rate folded with the WIMP velocity distribution is given by: 
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dir \ _ /^\ _ P(0) rn 
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where 
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(10) 



The event rate for the coherent WIMP-nucleus elastic scattering is given by 2^, 27, 13, 37 1 

p(0) TO 



R 



with 



fcoh{A, ^^r{A)) = 



lOOGcV 



A^r(A) 
Mr(p) 



A tcoh (1 + hcohcosa) 



(11) 



(12) 
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(13) 



with ^0 and a^^^ the scalar and spin proton cross sections Cspin the nuclear spin ME. In this 
work we will not be concerned with the spin cross section. 



The number of events in time t due to the scalar interaction, which leads to coherence, is: 



i?~ 1.60 10 



-3 



P(0) 



m 



ly 0.3Gel/cm-3 lKg280kms-i 10-6 pb 



fcoh{A,flr{A)) 



(14) 



In the above expression m is the target mass, A is the number of nucleons in the nucleus and (w^) 
is the average value of the square of the WIMP velocity (with n = 1). The quantity of interest to 
us is the quantity r = tcoh (1 + hcohcosa) , which contains all the information regarding the WIMP 
velocity distribution and the structure of the nucleus. It also depends on the reduced mass of the 
system. It is not difficult to show that: 



dr 
du 



where 



^{x,a)=l ydy J d^ J d(t)f{x,y,^,(j>,n) 



(15) 



(16) 



where F(u) is the nuclear form factor. The quantity ^a{y,a) enters since in some region of the 
velocity space the upper value of ^ is restricted so that the condition 



+ 2y^(l + 5 cos a) + 2(1 + 2(5 cos a + 46^) ^ 



nye 



is satisfied. 

The integration can be performed to yield: 



yesc rioiv.a) 

ydy J d^fi{x,y,£„n) 



with 



(17) 



(18) 
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(19) 



where Iq is the well known modified Bessel function. 

By performing a Fourier analysis of the function ^'(a;,a), which is a periodic function of a, and 
keeping the dominant terms we obtain the two amplitudes ^oio.\/u) and Hia^/u). Thus: 



dr dt dh 
du du du 



where 
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(20) 



dt 
du 
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(21) 



The total (time averaged) rate is given by: 



icoh — 



dtcoh 

du 



du 



(22) 



where 



( TIV )^ 

<^ detector threshold , Umax = -|- — ^ maximum WIMP velocity 



r = tcoh i'^ + hcoh cos a) , hcoh ^ -r— f ^^^^rfw (23) 



By including both ^!o(a^/u) and H{a^/u) we can cast the rate in the form: 

icoh J U^in 

The direction e of the recoiling nucleus is observed. 

In this case the directional differential rate is given by: 

(^) = '-^i^V<y) I ^^^^^fi^^.r^EZ-^^Si-^ - i.e)d^. (24) 
\ au / Amj^ J \/ (v ) "U ztt /irOw 

The factor of l/27r appears, since we are using the same cross section as in the non directional case, 
even though no angular integration is now required. The above coordinate system, properly taking 
into account the motion of the sun and the geometry of the galaxy, is not the most convenient for 
performing the needed integrations in the case of the directional expressions. For this purpose we go 
to another coordinate system in which the polar axis, Z, is in the direction of observation (direction 
of the recoiling nucleus) via the transformation: 

cos Q cos $ cos 9 sin $ — sin 8 

—sin^ cos*!' 
sin O cos $ sin Q sin <& cos Q 

In this coordinate system the orientation parameters and $ appear explicitly in the distribution 
function. In fact the numerator of the exponent of the M-B distribution (|10p becomes: 
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(5^+2y^ cos $ sin a sin 8(5— 2j/^ cos a cos 7 sin $ sin 8(5+2 cos a sin ^5—2y\/l — cos a cos <E> C0S7 sin((/))(5+ 
—2y\/l — sin a sin $ sin (j)6 + ~ 2?/ -y^l — cos (f) sin 8 [5 cos a sin 7 + 1)+ 

+2j/cos8 ^5 cos sin 7,^ + C + — cos $ cos 4) sin a — 5\/\ — cos a cos 7 cos((/)) sin<I>^ + 1 

This is clearly quite messy, but the constraint in the integration variables imposed by the Heaviside 
function becomes quite simple, namely the velocity in polar coordinates is specified by the angles 
< 6* < 7r/2 , < < 27r. 

The 5 function ensures that in the directional case the variables m, v and ^ obey the required 
relation. In our numerical calculation we found it more convenient to use y = v/vq and to express 
^ in terms of the other two, namely u and y. So one is left with two integrations, over (j) and y. To 
make the calculations tractable we made an expansion to first order in (5 = — = 0.135. Thus the 
distribution takes the form: 



^ 1 / y^ — y^ — cos (h sin 6 + 2a: cos 8 + 1 
][v,VE,Q,'^) = -^-thExp ^- 



[rm? + +2(5 cos q;x cos 7 sin <& sin 8 + y\/y^ ~ cos 4> sin 8 — ij sin 7 
+25cos(3!COs8 {^\/ y^ — cos 7 cos (f) sin $ — a::sin(7)^ b\J y"^ — x'^ sin a sin $ sin <^ 
—25 cos $ ^ — x'^ cos 8 cos (\) sin a + x sin 8 sin (a? — \/ — x'^ cos a cos 7 sin ] (25) 
The integration over the angle </> can now easily be accomplished yielding: 



• Time independent part: 



fo ix,y,Q) = 3 ^ Exp /o ^ (26 

Note that this part is independent of $. 

The modulation amphtude proportional to cos a 

/*'^x,2/,e,$ = ^^(5£;a;p -i^ ^ 27 



/ 2a/w2 — sin 9 \ 
[/o — ^- 2 '''^^ $ sin 6 — (x cos 9 + 1) sin 7) 



Hi — 2j i 



sin 9 



(cos 9 cos 7 sin <I> + sin 9 sin 7)] 

• The modulation amplitude which is proportional to sin a: 

/^(x,.,e,$) . -^^Exv (v1A:^^1^^^\ ,,3^ (28) 



I V2/ - X h I — 1 cos 9 + x/o I — 1 sm9 1 

Such an amplitude does not appear to leading order in non directional rate. Here it may 
become important near $ = or tt. 

In the above expressions Iq and Ii are the modified Bessel functions. 

From now on the needed integrations over y and u can be done numerically. Thus one obtains: 

vl/f (x,9,$)= / (a;,y,9,$)dy, z = 0,c,s; (29) 

H„,{{x, 9, $)) = (*f -(x, 9, $))2 + (*f (a;, 9, $))2 (30) 

The last function maybe used in obtaining the magnitude of the modulation amplitude. From now 
on we proceed as in the previous section, except that the obtained results are functions of 9 and $. 



^ V^'^'^"^**''^"^^' ^ V ^'^'^'^"^ (aV^)cosa + *f (aV^)sina] 



(31) 
or equivalently 

(^) = \/f«'^'N^™(«V^)cos(a + i9™7r) (32) 

where 'dm is the shift in the phase of the modulation (in units of tt) relative to the phase of the 
Earth, namely: 

= arctan V 33 

The event rate is still given by Eq. (HH) except that now: 
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^coh 



and 
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du 



(35) 



(36) 



K is a measure of the reduction in the event rate in directional experiments over and above the 
geometric factor of l/(27r). Furthermore 



he cos a -\- hs sin a 



1 



du 



(37) 



dir 



hm{coh) — \Jh?e + ^s, OLrn = arctan 



(38) 



It is clear that the quantifies and in Eq. ([57)1 are obtained from and ^'^ respectively. 
Clearly the quantities k, /im and a™ depend on the direction of observation. The range of the 
relevant angles needed to specify the line of recoil can be chosen to be : 

0<e<7r, 0<$<27r 

Sometimes in directional experiments one is interested in the angular distribution of events around 
the direction of observation, since the direction of the track may not be known precisely 
[ssj ]. To get an expression for this, we proceed as above by eliminating the variable y in terms of u 
and ^. Then for a given ^ we integrate over the energy transfer u. In this instance we will ignore the 
dependence of the event rate on (5, i.e. we will neglect the modulation. Thus we define the angular 
distribution of the expected events in time as follows: 



dRdir,coh 



1.60 10" 
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Ignoring the modulation effect we get: 
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(40) 



Partly directional experiments. 

In this case one can specify the line the nucleus is recoiling but not the sense of direction on it. 
The results in this case are obtained from those of the previous section via the replacements of the 
functions defined above by: 



ftix, y, 9) = ft{x, y, e) + /o*'-(x, y, TT - 9) 



(41) 



ff'ix, y, 9, $) = ff^ix, y, 9, $) + /^^(a;, y, ^ - 9, $ + tt), * = c, s 
The range of the angles now is: 

< 9 < tt/2, < $ < vr 



(42) 



(a) (b) 

FIG. 1: The function ^g"^ as a function of x and Q (y-axis in radians), (a) n = 1 and (b) n = 2. 



Some general results 

Before proceeding to specific applications involving specific nuclear targets it is instructive to make 
general observations, ignoring the fact that eventually x = a^/u, with a depending, among other 
things, on the specific target. First we note that the time average rate depends only on the polar 
angle Q (see Fig. [T]). From Fig. [T]we see that for large x the function ^'q"' peaks in the direction 
opposite to the sun's direction of motion , as expected for a M-B distribution. This, however is not 
true for small x. This is understood by observing that x appears explicitly as the coefficient of cos O 
in the exponential of Eq. (|26l ) The situation for the modulated amplitude H„i is more complicated 
since it also depends on the angle Therefore we will not discuss this case here. We only mention 
that the modulation here in the case of the n = 2 is suppressed relative to that for the standard 
n = 1 case, independently of the angle of observation. This is expected in view of the results for the 
standard (non directional) case [36|. 

SOME APPLICATIONS 

In this section we are going to apply the formalism of the previous section in the case of two 
popular targets: i) The hght target appearing in CS2 involved in DRIFT [20*1 and ii) The ^^''l 
target, which has been employed in the DAMA experiment [39, 40]. We will not consider energy 
thresh hold effects and we will ignore quenching factor effects. We will consider only the coherent 
mode, but the results obtained for the functions t, k, hm and am are not expected to be radically 
modified, if one considers the spin mode. 

The light target CS2 

The nuclear form factor employed was obtained in the shell model description of the target and 
is shown in Fig. [21 The parameter t entering the non directional case is shown in Fig. [3] 




(a) (b) 



FIG. 2: (a) The form factor F^{u) for ^^S employed in our calculation with u = Q /Qo, Q the energy transfer 
to the nucleus and Qo = 404 keV. (b) The same quantity a function of the energy transfer Q. 




FIG. 3: The quantity tcoh is shown for Qmin ~ with (a) n — 1 and (b) n = 2. 



We begin our analysis of the directional signal by computing the differential rate with respect to 



^, which is used in simulating the experiments [2i 
and [5] 



\m 33 



Our results are presented in Fig. S] 
With the above angular distribution we can obtain the average value of ^, < ^ >. Thus we 
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FIG. 4: The differential rate dr/d^ for ^^^S in as function of ^, the cosine of the angle between the line 
of observation and the line of recoil. The thick solid, fine solid, short, and long dash correspond to = 
7r/4, 7r/2, 37r/4 and vr respectively, (a) directional in a given sense for n = 1 in the case of a WIMP mass of 
10 GeV. (b) The same as in (a) for n = 2. (c) the same as in (a) for 30 GeV. (d) the same as in (b) for 30 
GeV. 



get the results shown in Fig. [51 We see that the average value of ^ does not change vary much in 
going from the rt = 1 to the n — 2 case. We clearly see that this value is much higher in the case 
the observation is made opposite to the sun's direction of motion (0 = tt). 

We continue our analysis by calculating the parameter Hcoh discussed above. Our results are 
shown in Figs [TlfTOl From these figures we see that in the case of directional experiments the 
maximum value of k is attained when the observation is made opposite to the sun's direction of 
motion. This is expected in the case of the M-B distribution. We have seen above that < ^ > also 




(c) (d) 



FIG. 5: The same as in Fig. [l] (a) directional in a given sense for n = 1 for 100 GeV. (b) The same as in 
(a) for n — 2. (c) the same as in (a) for 250 GeV. (d) the same as in (b) for 250 GeV. 




FIG. 6: The quantity < ^ > is shown for Qmin = with (a) n = 1 and (b) n = 2. Otherwise the notation 
is that of Fig. g] 



attains its maximum value. When the sense of direction is not observed the maximum event rate is 
again attained when the observation is along the line of the sun's motion. 

We will next discus the modulation effect, i.e. the parameters hm, shown in Figs 12511^ and am 
shown in Figs fTMTSl We clearly see that the results depend on the WIMP mass. From the figures 
rmfHl we see that the modulation for n = 2 is quite small, reminiscent of a similar result in the case of 
the non directional case ^3^. Since the exhibited pattern is otherwise similar to the standard n = 1 
case we will limit our discussion to the n = 1 case. The modulation amplitude, h^, ~ 0.1, in the 
most favored direction O = tt, i.e. opposite to the sun's direction of motion, is small, but still much 
larger than that encountered in the non directional case. The absolute maximum hm ~ 0.4 seems 
to occur when $ = tt, i.e. inwards towards the center of the galaxy, at an angle O — 0.5 ~ 7r/6. 
A smaller value of 0.3 is attained when one is looking radially outwards from the center of the 




(c) (d) 
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FIG. 7: The parameter k defined in tlie text for a iigiit target in tiie case of a WIMP mass of 10 GeV. (a) 
Directional in a given sense for n = 1. (b) The sum of both senses for n = 1. (c) The same as in (a) for 
n — 2. (d) The same as in (b) for n — 2. 




(c) (d) 
> radians 



FIG. 8: The same as in Fig. [3 for a WIMP mass of 30 GeV. 



galaxy at Q — 1.0 ~ 7r/3. Another maximum /i,„ « 0.2 is attained, when 0«7r/2, i.e. in a plane 
perpendicular to the sun's direction of motion, and $ = 7r/2, along the line perpendicular to the 
galactic plane. 

When the sense of direction is not observed the modulation picture is quite simple. One sees that 
the modulation is quite small, when the line of observation is along the sun's motion and reaches a 
maximum value, hm ~ 0.3 on the plane perpendicular to the sun's velocity regardless of 

Regarding the phase am we see that, as expected, it is zero at = or tt, reminiscent of the 
non directional case. For <f> — Tr/2 or 37r/2 this phase is zero for all Q. This means that in all these 
cases one has the standard behavior of the modulation (maximum around June 2nd). For <I> = 0, 
outwards from the center of the galaxy, the phase decreases from zero to n/2 depending so long 
as O < 0. During this period the maximum precedes the phase of the Earth. When O > 7r/2, it 
changes sign (the maximum drags behind the phase of the Earth). In other words the maximum 
occurs in the spring for 6 > •7r/2 and in the autumn for Q < 7t/2. The opposite is true when $ = tt. 




(c) (d) 

— > radians 



FIG. 9: The same as in Fig. [7] for a WIMP mass of 100 GeV. 




(c) (d) 

O ^ radians 



FIG. 10: The same as in Fig. [7] for a WIMP mass of 250 GeV. 



In the case when both senses of WIMP recoil are considered the modulation drags behind the 
phase of the Earth, a maximum in the spring, for <i> = 0, tt (radially in the galaxy). It attains a 
maximum in the fall for = p/2, 37r/2 (perpendicular to the plane of the galaxy) 

Anyway the above complicated pattern of the shift in the location of the maximum, which depends 
on the LSP mass, is not expected to cause problems in the analysis of the experiments. The periodic 
nature of the phenomenon is sufBcient. 



The intermediate mass target 127I 

The nuclear form factor employed was obtained in the shell model description of the target and is 
shown in Fig. [111 We first exhibit the previously obtained [3^ parameter tcoh in Fig. We have 
a suppression as the WIMP mass increases (partly counteracted by the factor /i^ / mwiM p entering 



(a) (b) 
G > radians 




FIG. 11: The parameter hm, defined in the text, as function of the polar angle Q for a light target in the 
case of a WIMP mass of 10 GeV. The fine solid, thick solid, dashed and long dashed curves correspond to 
"1> = 0, 7r/2,7r, and 3n/2 respectively (the cases $ — 7r/2 and 37r/2 cannot be distinguished). Note that, 
when the sense is not distinguished, the maximum angle of the line of observation with the sun's direction 
of motion is n/2, while the maximum azimuthal angle is vr. Now the thick solid, fine solid, dashed and long 
dashed curve correspond to "I> = 0, 7r/4, 7r/2 and 37r/4 respectively (the cases $ = 7r/4 and 37r/4 cannot be 
distinguished), (a) The directional rate in a given sense for n — 1. (b) The sum of both senses is shown for 
n = 1. (c) The same as in (a) for n — 2. (d) the same as in (b) for n = 2. 




FIG. 12: The same as in Fig. [TT]for a WIMP mass of 30 GeV. 



the event rate and not included in tcoh)- The behavior of the parameter k is exhibited in Figs [2Tll24l 
One can see that for relatively light WIMP in the directional case k has a similar behavior with that 
encountered above for a light target. It peaks in the direction opposite to the sun's velocity, i.e at 
O = TT, as expected. For heavy WIMP masses the maximum value is a factor of 2 smaller than that 
for a light target and it occurs at lower polar angles 8. When both senses are considered for n = 1 k 
peaks along the sun's line of motion (by definition here = 0). The situation is reversed for heavy 
WIMP, i.e. it attains a maximum at O = tt/2. This reversed pattern holds for n — 2, except for 
very small WIMP masses. Thus, if one could independently establish that the velocity distribution 




FIG. 13: The same as in Fig. [TT]for a WIMP mass of 100 GeV. 




FIG. 14: The same as in Fig. [TT]for a WIMP mass of 250 GeV. 



is of M-B type, this signature can be used in constraining the WIMP mass. 

We will next discus the modulation effect, i.e. the parameters ft,,,,,, shown in Figs [25ll28l We 
clearly see that the results depend on the WIMP mass. From the figures [25ll28l we see that in the 
directional case the modulation is similar to that for a light target for small WIMP masses. This 
is not surprising, since, then, the reduced mass the same for both targets. For heavy WIMPs the 
modulation goes through a small value at around Q — {2/3)tt regardless of the angle $. The 
maximum value of hm is also a bit smaller. 

If both senses are counted the modulation amplitude becomes essentially independent of $ and 
attains the maximum value of hm = 0.1 — 0.3 depending on the WIMP mass. Note that in this case 
for heavy WIMPs one has a secondary maximum and a minimum at ~ (3/4) w 7r/4, .i.e. when the 
recoil occurs in the plane perpendicular to the sun's velocity. Again the absolute maximum occurs 
at = vr, i.e. along the sun's motion. 

We will not show the phase am , since the obtained pattern is similar to that encountered for a light 
target. Before concluding this section we should mention again that in the case of ^^^I one may 
have a contribution due to the spin cross section. As we have already mentioned, the controversy 




FIG. 15: The parameter a™ (in units of tt), defined in the text, as function of the polar angle Q for a light 
target. Otherwise the notation is the same as in Fig. 1111 




regarding the DAMA experiment, employing Nal target, may be attributed to this interaction, 
which does not enter in experiments involving even nuclear targets. This possibility is currently 
under study, including our own realistic (static) spin ME and spin form factors. The quantities 
t, h, dr /d^, K and hm do not depend on the spin ME, they only depend on the adopted spin form 
factors. We do not, however, expect the factors dr/d^, k and /i„i, which are the main subject of this 
work, to be radically different from those presented here for the coherent cross section. 



CONCLUSIONS 



We have seen that, given a sufficient number of events, the directional experiments, in which 
one measures the direction of the nuclear recoil, provide an excellent signature to discriminate 
against background. Some of these features persist even if the sense of motion of recoils along their 
line of motion cannot be measured. The predictions depend, of course, on the assumed velocity 




FIG. 17: The same as in Fig. [15] for a WIMP mass of 100 GeV. 




FIG. 18: The same as in Fig. [15] for a WIMP mass of 250 GeV. 




FIG. 19: (a) The form factor F^{u) for ^^'^I employed in our calculation with u = Q/Qo, Q the energy 
transfer to the nucleus and Qo = 64 keV. (b) The same quantity a function of the energy transfer Q. 



distribution. In the present work we selected to work with a M-B distribution: (i) the traditional 
one with characteristic velocity that of sun around the center of the galaxy and (ii) a variant obtained 
when dark matter is coupled to dark energy via a scalar field yielding an increase in the gravitational 
field for dark matter [36|. In the latter case the characteristic velocity in the WIMP distribution 




increases by a factor or n > 1. 




In the most favored direction, opposite to the sun's direction of motion, the event rate is ~ ^ 
down from the standard non directional experiments. The modulation amplitude in this direction 
hm depends on the WIMP mass. For a light target it ranges between 0.05 and and 0.1 depending 
on the WIMP mass. For a heavy target it is somewhat reduced, hm — 0.02 — 0.07, but it remains 
higher than that expected in the non directional experiments. It is also characterized by a definite 
sign (maximum around June 3nd). Higher values h^n — 0.25 — 0.35, yielding 50 — 70% difference 
between the maximum and the minimum rates, are expected by a judicious choice of the direction 
of observation. Thus quite large asymmetries with seasonal dependence are predicted. The time of 
the maximum is also direction dependent, so it cannot be mimicked by seasonal variations of the 
background. 

In partly directional experiments, i.e. experiments in which the sense of motion of the recoiling 
nucleus is not determined, one no longer can measure asymmetries. We still find, however, that the 
expected event rate is maximum along the line of motion. Equally large modulation amplitudes hm 
are predicted with a seasonal variation, which, again, cannot be mimicked by seasonal variations in 
the background. 





The modulation amplitude is decreased, if dark matter is coupled with dark energy, in a fashion 
analogous to the non directional case. This is not unexpected, since the ratio of the velocity of the 
earth to the characteristic velocity becomes smaller. This suggests that one should test whether the 
above conclusions hold, by considering other velocity distributions. Thus one may be able to infer 
the velocity distribution, when the experimental data become available. 

The directional experiments are quite hard. It is encouraging that the angular distribution is 
such that the average angle is given by < ^ >=< cos9 >— 0.5 — 0.8 depending on the angle 
0, the polar angle of the direction of observation from the sun's direction of motion. Anyway in 
experiments planned by DRIFT [i^l one can register events recoiling in all directions. If some 
interesting events are found, one may further analyze them by the direction of recoil to discriminate 
against background. 



(a) (b) 
G > radians 




FIG. 27: The same as in Fig. [25] for a WIMP mass of 100 GeV. 




FIG. 28: The same as in Fig. iSlfor a WIMP mass of 250 GeV. 
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